FLAVOR CHANGING KAON DECAYS FROM HYPERCP: 
MEASUREMENTS OF THE K± -» it^n+yT BRANCHING RATIOS 



E. Craig Duke^ 
Physics Department, University of Virginia, 
Charlottesville, VA 22901, U.S.A. 



The Fermilab HyperCP collaboration is making precision studies of charged hyperon and 
kaon decays, as well as searches for rare and forbidden hyperon and kaon decays. We report 
here on measurements of the branching ratios of the flavor-changing neutral-current decays: 
K — > n ± u + , and compare our results to theoretical predictions. This is the first observa- 
tion of the K~ — » 7T~ u + fi~ decay. 



1 Introduction 



The radiative, flavor-changing neutral-current decays of the type — ► 7r ± Z + / _ were originally of 
considerable interest because, as tree-level forbidden decays in the standard model, they were im- 
portant systems in which to search for new physics. The interest in these decays has not waned, 
despite the lack of any indication of new physics, but the focus has shifted to using them as a test- 
bed for techniques in calculating quantum-chromodynamic (QCD) corrections to weak decays. 
Expected to proceed via one-loop level electroweak box and penguin diagrams, long-distance 
effects are thought to play a dominant role in these decays. Hence chiral perturbation theory 
(ChPT) has been employed to calculate both the rates and the form factors of these decays. 
An early ChPT calculation which parameterized the rates and form factors in terms of a single 
parameter, has recently been supplemented by a higher-order, 0(p e ), ChPT calculation 

"Representing the HyperCP collaboration: A. Chan, Y.C. Chen, C. Ho, P.K. Teng, (Academia Sinica); 
W.S. Choong, G. Gidal, Y. Fu, P. Gu, T. Jones, K.B. Luk, B. Turko, P. Zyla, (Berkeley and LBNL); C. James, 
J. Volk, (FNAL); J. Felix, (Guanajuato); R.A. Burnstein, A. Chakravorty, D.M. Kaplan, L.M. Lederman, W. Lue- 
bke, D. Rajaram, H.A. Rubin, N. Solomey, Y. Torun, C.G. White, S.L. White, (IIT, Chicago); N. Leros, J.P. Per- 
roud, (Lausanne); H.R. Gustafson, M.J. Longo, F. Lopez, H.K. Park, (Michigan); K. Clark, M. Jenkins, (S. 
Alabama); E.C. Dukes, C. Durandet, R. Godang, T. Holmstrom, M. Huang, L. Lu, K.S. Nelson, (Virginia). 



which parameterizes them in terms of two variables: a+ and &+□ This last calculation makes a 
model-independent prediction that the ratio 1Z = T(K + — > n + fi + fi~) /T(K + — > ir + e + e~) should 
be greater than 0.23, a result of rather broad validity since it assumes only that the electron and 
muon channels are described by the same electromagnetic form factor and that this form factor 
does not have singular behavior. A previous calculation, using a different technique predicts a 
similar result! 

Although first seen in 1975, with the observation of some forty K + — > ■K + e + e~ (K^ ee ) events^ 
only in the last several years have precision measurements been made of the K^ ee branching ratio 
and form factor in two experiments at the Brookhaven National Laboratory (BNL): BNL-777^ 
and BNL-8650 The most recent of these measurements comes from BNL-865. Using a sample of 
10 300 events they have firmly established the vector nature of the decay, measured the slope of 
the form factor to 9%, and the branching ratio to 5%. The K + — > 7r + fi + fi~ (^Jlm) branching 
ratio measurement has been more difficult to make. First observed in 1997 by BNL-787^1 they 
found a branching ratio about 2a less than that expected, given the prediction of 7Z mentioned 
above and the BNL-E865 measurement of the K+ ee branching ratio. The K^L^ branching ratio 
has recently been measured by BNL-865 with a sample of some 400 eventsta This latter result 
agrees with the theoretical predictions and disagrees with the BNL- 787 result by 3.3<r. Hence it 
is vital that this discrepancy be resolved and to determine if indeed the theorists predictions of 
1Z are borne out. 



2 The HyperCP Apparatus 

The HyperCP experiment was designed and built primarily to investigate CP violation in 
H 1 * 1 — ► Air^ hyperon decays. A plan view of the HyperCP spectrometer is shown in Fig. Q. A 
charged-secondary beam was produced by steering an 800-GeV/c proton beam onto a 2x2 mm 2 - 
Cu target. A curved collimator with a 4.88 ^sr solid angle acceptance and embedded in a 6-m 
long dipole magnet deflected charged particles, with an average momentum of about 170 GeV/c, 
up at a mean angle of 19.5 mrad. Following a vacuum decay region was a conventional magnetic 
spectrometer employing high-rate, narrow-pitch wire chambers. Two hodoscopes on either side 
of the channeled beam eminating from the collimator exit, called the Left-side and Right-side 
hodosopes, were the basis of most of the physics triggers. A hadronic calorimeter was used 
for triggering on the proton from A — > pir decays, and was not used in the analysis reported 
here. At the rear of the spectrometer was a muon system consisting of two identical detectors 
positioned on either side of the channeled beam. The detectors had three proportional-tube 
stations, each behind 0.81 m of iron absorber, providing x and y measurements. At the rear of 
the muon detectors were two hodoscopes, with vertical and horizontal counters, used to trigger 
on events. 

The polarities of both the Hyperon and Spectrometer magnets were periodically flipped to 
change the sign of the secondary beam; roughly two-thirds of the running time was spent on 
the positive polarity. A typical primary beam intensity of 6xl0 9 protons per second gave a 
secondary beam rate of 13 MHz at the exit of the collimator. Data were taken in two runs, one 
in 1997, another in 1999, with a total of 231 billion events recorded on magnetic tape, of which 
555 million were fully-reconstructed — > it zil, K + 7t~ (Kjj^) events. The analysis reported here 
is from the smaller 1997 data set containing 58 billion events. 



3 Analysis 

Rather than measuring absolute branching ratios, which would have required the somewhat 
difficult task of counting the number of kaons entering the apparatus, we measured the 
branching ratios relative to the well-measured (and far larger) branching ratios. We call 




Figure 1: Plan view of the HyperCP spectrometer. 



Kff^ the normalization modes and K^„„ the signal modes. Since the topologies of both the 
signal and normalization modes were quite similar, measuring the branching ratios by this means 
minimized systematic errors. 

The triggers used for both signal and normalization modes were both quite simple. For the 
normalization modes the triggers required a coincidence of signals in both the Left-side and 
Right-side hodoscopes. The triggers for the signal modes required this same coincidence along 
with a coincidence of charged particles in the Left-side and Right-side muon hodoscopes. 

The branching ratios for were determined using the following formula: 



1 N obs A f tri 9 f sel R(K ± — > tt^tt+tt- 

[ 77 V" V J_ 900iV o6s A trig sel p+p- ' " 



where: N obs are the numbers of observed events, A are the geometric acceptances of the triggers, 
e sel are the event-selection efficiencies, and e tri9 are the trigger efficiencies; for signal (K^ ) and 
normalization (K^ nn ) modes. The number 200 is the prescale factor for the normalization trigger 
and e M+/x is the muon identification efficiency. The values of the parameters used in Eq. [I] are 
given in Table [l] below. Note the generally large values of the efficiencies and acceptances. The 
following discussion explains how these values were obtained. 



Table 1: Parameters used to calculate the branching ratios. 
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Data for both signal and normalization modes were selected in a similar manner. A three- 
track topology — two on the left side of the spectrometer and one on the right — was required 
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Figure 2: The 7r/x/i invariant mass for positive and negative beam data after cuts. The large peak below the K 
mass is from K^ nn events with either in-flight decays of the pions or pion punchthrough. The insets show the 
signal region used to establish the number of events. 



in both cases, with the signal tracks having the additional requirement of a muon track in both 
the Left- and Right-side muon detectors. A muon track was defined as a track stub in a muon 
detector which included at least two of three hits in both x and y views in the muon proportional 
tubes with corresponding in-time hits in the muon hodoscopes. Two of the three pion tracks 
for the normalization modes were required to point to the fiducial areas of the muon detectors. 
A good fit to a common vertex was required of the three tracks. This was done by first fitting 
the tracks to a common vertex and requiring that the average transverse separation of the three 
tracks at that point be less than 2 mm and the x 2 /d.o.f of the vertex-constrained fit be less than 
2.5. This cut out most of the background from hyperon decays and uncorrelated tracks. Other 
rather loose cuts were applied to the parent: on its momentum, target point-back, and decay 
vertex. The invariant mass of the signal and normalization modes after these cuts is shown in 
Fig. ||[ The number of signal events, as determined by maximum-likelihood fits in the region 
470 MeV/c 2 < < 520 MeV/c 2 , is given in Table |l|. The invariant mass of normalization 

data sample is shown in Fig. ||. 

From the combined positive and negative data samples we measured the relative trigger 
efficiency, e^^/e^"-?, and the dimuon selection efficiency: 6^+^-, given in Table Q 

From Monte Carlo we determined the acceptances and the event selection efficiencies for 
signal and normalization modes. The Particle Data Group (PDG) parameters describing the 
decays were usedu The Kjjjr„„ decay parameters are not as well known. In the Monte 
Carlo simulation we assumed a vector interaction, as has been demonstrated in K+ ee decays 
by the BNL-865 experiment^ and is consistent, although not nearly as well established, with 
their measurement! For the slope of the form factor which describes the dilepton invariant 

mass spectrum, Mu, for the K^„„ decays we averaged the values from the BNL-777 K+ ee and 
BNL-865 K+ ee and K£ measurements to obtain 5 = 2.08 ±0.12: the form factor being given by 
4>(M^) = 0o (1 + dM^/m 2 K+ ), where 5 and 4>o are a slope parameter and a constant, respectively, 
and m K + is the mass of the charged kaon. It is important to note, as can be seen in Fig. |||, 
that our acceptance in the dimuon invariant mass was quite uniform and hence we were not 
sensitive to even large deviations from the average value of the slope parameter 5. The geometric 
acceptances are given in Table The event selection efficiencies, also measured using Monte 
Carlo, are quite high, and are also given in Table |l|. 

A word on errors. Our measurements were dominated by the statistical uncertainties in 
determining the number of Kijjr events: giving a 10% error for the combined branching ratio 




Figure 3: The -k~-k + tt~ invariant mass. 
The dashed and dotted curves, bi and D2, 
are respectively linear and quadratic fits to 
the background in the side-bands. 



Figure 4: The dimuon invariant mass for (a) positive- and (b) 
negative-beam data after cuts. The solid lines are the Monte 
Carlo predictions using the measured values of the slope pa- 
rameters. The dashed and dotted lines represent Monte Carlo 
predictions with ±3cr variations, respectively, from the measured 
value of the slope parameter. The curves show the acceptance. 



measurement. Sources of systematic error were carefully studied, including: targetting vari- 
ations, magnetic field variations, the difference between data and Monte Carlo, uncertainties 
in the parameters describing the K^„„ and decays, and uncertainties in the trigger and 

muon identification efficiencies. The total systematic errors were found to be less than half the 
statistical errors. They were dominated by two sources roughly equal in magnitude: the error in 
the background under the signal peak (3.6%) and the error in the determination of the relative 
trigger efficiencies (3.1%). 



4 Results 

What do we find for the branching ratios? Using Eq. [l] we get: B(K+ ) = (9.7=tT.2±0.4) x 1CT 8 
and B(K~ ' ) = (10.0 ± 1.9 ± 0.7) x 10~ 8 , where the first and second errors are statistical 
and systematic, respectively. The two results are in agreement with each other, indicating no 
evidence of CP violation: 
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No surprise here: theoretical expectations for the CP asymmetry are orders of magnitude 
smallerli^ Assuming CP symmetry is good, we can combine the two measurements to give: 

B(K^) = [9.8±1.0(stat)±0.5(syst)]xl0 -8 . 

This result is consistent with the BNL-865 measurement, but 3.2a higher than the BNL-787 
result (see Fig. ||). It also, when combined with the PDG value for B(K^ ee ), is consistent with 
the theoretical predications of 1Z = T (K^^) /T (K^ee) ■ 

The HyperCP collaboration is analyzing the larger 1999 data set, which should give about 
a factor of four increase in statistics, and which has quite different systematics: a different 
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Figure 5: Experimental measurements of the 
branching ratios. The two BNL points refer to mea- 
surements of B(Knuu)- 
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Figure 6: The nfifj, invariant mass from the 1999 
HyperCP data set. 



muon detector, Monte Carlo, and tracking program. Figure ^ shows the 7Tfifi mass peak from a 
preliminary analysis of the combined K£ and K~ 1999 data sets. By the end of the decade, 
new experiments such as CKM at Fermilab, should produce data sets for both the K^ ee and 
modes that are on the order of 100,000 events, which will allow precision tests of ChPT 
predictions to be made. 

We finally note that these results have been recently published in Physical Review LettersS 
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